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ABSTRACT 

We point out an important source of error in measurements of extragalactic Ha emis- 
sion and suggest ways to reduce it. 

The Ha line, used for estimating star formation rates, is commonly measured by 
imaging in a narrow band and a wide band, both which include the line. The image 
analysis relies on the accurate removal of the underlying continuum. We discuss in 
detail the derivation of the emission line's equivalent width and flux for extragalactic 
extended sources, and the required photometric calibrations. We describe commonly 
used continuum-subtraction procedures, and discuss the uncertainties that they intro- 
duce. 

Specifically, we analyse errors introduced by colour effects. We show that the 
errors in the measured Ha equivalent width induced by colour effects can lead to 
underestimates as large as 40% and overestimates as large as 10%, depending on the 
underlying galaxy's stellar population and the continuum-subtraction procedure used. 
We also show that these errors may lead to biases in results of surveys, and to the 
underestimation of the cosmic star formation rate at low redshifts (the low z points 
in the Madau plot). We suggest a method to significantly reduce these errors using a 
single colour measurement. 
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1 INTRODUCTION 

Physical properties of galaxies can often be derived from 
spectroscopic observations, which allow the detection and 
characterisation of various emission lines. However, the com- 
plexities related with observation and reduction procedures, 
make spectroscopic surveys considerably more difficult to 
pursue than imaging surveys using narrow-band filters. 
While such imaging surveys detect only lines whose red- 
shifted wavelength fall inside the filter's band, they have the 
advantage of detecting spatially extended line emission down 
to low surface brightness levels. This imaging technique is, 
therefore, useful for evaluating the full spatial extent and dis- 
tribution of the line-emission in nearby galaxies (Kennicutt 
et al. 2008), as well as to search and identify high-redshift 
line-emission galaxies in studies where many objects are tar- 
geted (e.g., Thompson, Mannucci & Beckwith 1996; Nilsson 
et al. 2007). 

The Ha line provides a large part of our knowledge 
of star formation (SF) activity in galaxies and is the fo- 
cus of this paper. Observing a galaxy with a filter centred 
at rest-frame Ha can, with reasonable assumptions, provide 
information on the amount and spatial distribution of the 
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emitted ionising radiation using relatively modest integra- 
tion times on small telescopes, thus making large galaxy 
surveys practical. 

Since the flux measured by narrow-band imaging in- 
cludes both the spectral line and the continuum emission, 
it is necessary to measure the continuum contribution. In 
many cases this is done by imaging the target in a band 
significantly wider than the narrow band (e.g., Waller 1990; 
James et al. 2004; Kennicutt et al. 2008; Brosch, Spector & 
Zitrin 2011). A scaling factor between the continuum counts 
of the wide and narrow band images (Wide to Narrow Con- 
tinuum Ratio - WNCR) is then determined. From its value 
the line's counts can be extracted. 

This paper discusses the derivation of emission line 
properties from wide and narrow band imaging, focusing on 
sources of measurement errors. In Section [5] we derive the 
relations required to calculate the equivalent width (EW) 
and the flux. Section [3] discusses methods for obtaining a 
calibrated measure of the flux. In Section [3] we derive the 
effect of the WNCR estimation error on the measurement 
errors of the line's equivalent width and flux. In Section[5]we 
review some commonly used methods to measure WNCR. 
Then, in Section [6] we analyse the measurement errors that 
originate from colour effects when using these methods, and 



© 2011 RAS 



2 0. Spector, I. Finkelman and N. Brosch 



quantify these errors for the case of Ha measurements in the 
local Universe. 



2 DERIVATION OF EMISSION-LINE 
PROPERTIES 

Line emission from a galaxy or a region within a galaxy 
(target) can be measured by observing it through a narrow- 
band filter (N) and a wide-band filter (W), where the N 
band contains the emission-line, and the W band may or 
may not contain it. In this section we derive the method 
to calculate the emission-line's equivalent width (EW) and 
flux (Fune) from such observations, and discuss the different 
assumptions made in this derivation and their implications. 
The derivation and discussion are general, and can be ap- 
plied to Ha measurements as well as to other emission lines 
measured in a similar way. 

The response of the electro-optical system can be ex- 
pressed as the integral of the observed flux density over 
wavelength, weighted by the transmittance functions of the 
filter and of the atmosphere and by the response function of 
the rest of the electro-optical system: 

f oo 

cp Sjv = / Fa T atroiJV (A) T^ (A) R A dX 



where: 



cps w = FxT atm , w (X)T w (X)RxdX (f) 
'o 



cps^, cps w are the measured count rates of the narrow- 
band (N) and wide-band (W) filters (respectively) in instru- 
mental units - typically ADU s _1 (Analog to Digital Units 
per second); 

Fa is the spectral flux density of light from the observed 
target reaching the telescope in absence of an atmosphere 
(typically in erg s _1 cm -2 A ); 

Tjv(A), Tw{X) are the unitless transmittance functions of 
the N and W bands, respectively, ranging from (fully 
opaque) to 1 (fully transparent); 

T a t m ,jv(A), T a tm,w(X) are the atmospheric transmittance 
as function of wavelength, including effects of weather, ele- 
vation and airmass of observation, when the N and W bands 
(respectively) were imaged (unitless, ranging from to 14); 

Ra is the responsivity as function of wavelength of the rest 
of the electro-optical system (i.e. the telescope and sensors, 
excluding the transmittance effect of the filters) typically in 
ADU erg" 1 cm 2 ; 

A is the wavelength (typically in A). 

This paper considers the case where the measured line 
is the only significant spectral feature in the N and W bands, 
i.e. that the target's flux can be divided into line and con- 
tinuum (cont) emission. For the case of Ha, this means that 
some adjustments should be made, if the [Nil] lines flanking 
it carry significant flux. 

The count rate from the observed target can be ex- 
pressed as the sum of the line and continuum contributions: 



cps^ — cps Ncont + cps N j ine 

C P S W ~ C P S W,cont + C P S W,line 



(2) 



,cont = / F\ cont T a t m; iv(A) Tjv(A) Ra d\ 

p oo 

line / ^X,line T-atm,N (A) Tjv(A) Ra dX 

Fa, cont T a tm,w(X) Tvr(A) Ra dX 



cpsjv 

CPSjv 
C P S W,co 



C P 8 Wline = / F X ,UneT a tm,w(X)Tw(X)R\dX (3) 
JO 

where: 

cps Njine , cps N cont axe respectively the line and contin- 
uum contributions to cps^; 

c P s w,Une> c P s w,cont are respectively the line and contin- 
uum contributions to cps^; 

Fa, ;ine, Fa. cont are respectively the line and continuum 
contributions to Fa- 



Since the width of the emission-line is typically very 
narrow in comparison to both the N and W bands, its profile 
can be approximated as a delta function at the centre of the 
emission line: 



F \ } line ~ F^ ne <5(A A;j ne ) 



(4) 



where: 



Anne is the central wavelength of the emission line (typi- 
cally in A); 

Fime = j °° F\,Une dX - the flux (irradiance) of the emis- 
sion line (typically in erg s _1 cm 2 ). 

Using this approximation, the line contribution is sim- 
ply: 

C P S N,line — ^ line T ' a tm,N (Xune) Tjv(A; ine ) Ra(Ak„ 6 ) 
C P S W,line — Fiine T atm.W {Xline) Tw(Xline) RA(Aiine) (5) 

The goal of the measurement is to find the emission- 
line's flux, Fane, and equivalent width, EW. The line's 
equivalent width is defined as: 



EW 



F/ine 



Fa, cont (Azine) 

while, the flux can be extracted from ([5]): 

°P S JV,(ine 



Fi - — 

- 1 - Line — 



T a t m ,iV (Ajine) Tjv(AHne) RA(Ajine) 



(6) 



(7) 



The term T atm ,jv(Ai„ le ) ■ RA(A i4 ne) is obtained by pho- 
tometric calibrations and is discussed in section [3] The N- 
band's transmittance, Tjv(A;i ne ), measured directly or pro- 
vided by the filter manufacturer, should be corrected to com- 
pensate for effects of the telescope's converging beam, since 
the transmittance curves of interference filters (the type of 
filters usually used) vary with angle of incidence. 

We now derive the remaining terms: cpsjy Une and 
Fa, cont{Xune)- Using the two parts of ([5]), we can write 
cps w line in terms of cps N Une (based on the <5 function line 
approximation) : 

T a tm,w(Aime) Tw(Azjne) , , 

C P S W,line — C P S N,line Tp 7T TTp 7T T \°) 

i-at-m,N {Aline) *-N\Mine) 

Using this and ((2| we can express cps w cont as: 
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^ T^atm.W (Xline) Tw(Xline) 

C P S W,cont = C P S W ~ C P S N,line 77^ 7T \TZ 77 T 

± atm,N \s\line) -L N\Mine) 

(9) 

For the derivation of cps^ line we define the Wide to 
Narrow Continuum Ratio (WNCR) to be the ratio between 
the count rate contributed by the continuum in the W band 
and the count rate contributed by the continuum in the N 
band: 

cps w< 

CO li t 



WNCR : 



C P S N,COTl 



(10) 



In sections [4j [5] and [6] we describe the methods most 
commonly used to measure WNCR, and analyse their effect 
on the accuracy of EW and Fi me measurements. 

From |[2|, (JU and (flQ)| we obtain an expression for 

C P S N,line ] 



CPSjv.iine = CPSjv ~ CpS N cont = C PSjv 



cps v 



cps^ 



WNCR 

t aim, W inline) ^ W (Mine) 
C P S W ~ C P S N,line TP ' 77 TT 77 T 
' -L atm,N (Mine) ^N(Mine) 



WNCR 



Extracting cps Nline yields: 



(11) 



C P S N,line 



,cps N 



cps w \ 
' WNCR/ ' 

1 atm,W inline) Tw inline 



WNCR T a 



(12) 

Note that this expression includes the atmospheric 
transmittance at A; ine for both bands: T a tm,w(Xiine) and 
Tatm,N (Mine)- These can be estimated from a photomet- 
ric calibration in the N band, when available. If the W and 
N exposures were taken at very similar airmasses and at- 
mospheric conditions, these can be approximated as equal 

^ aim, W (Mine) ^ -| 

T Btm , _ 

Typically, the N band is chosen so that Tjv(Azine) is as 
high as possible (Xu n e is close to the N band's peak), and 
is not significantly smaller than Tw(Xune)- The W band is 
typically chosen to be significantly wider than the N band, 

therefore WNCR > 1, and also wnC&t'^a ) L ThiS 

means that the error due to small deviations from the above 
approximation will typically be minor. It will decrease the 
closer \n ne is to the N band's peak and the wider the W 
band is compared to the N band. 

In order to derive Fx, C ont(Xune) we approximate F x , CO nt 
to be constant within the N band. This is valid if the N 
filter is very narrow, and has negligible "out of band leaks" 
(a requirement which is not trivial for very narrow bands). 
If F\,cont is constant within the N band, it can be taken 
out of the integral for cps^ cont in ((3]) and be replaced by 

F X.cont inc. I ■ 



C P S N cont — cont{\line) / T atm , JV ( A) TjV ( A) Ra dX (13) 



Using ([2]), we get: 



t(A|. 



cps N - cps N li 



J °° T a t m> jv(A) Tjv(A) Ra dX 



(14) 



Replacing F \, c <mt{hme) and F Hne in © with the ex- 
pressions from (fT4)l . and using equation (0 yields: 



EW^ 



C P S N,line 



J °° T atm ,jv(A) Tjv(A) Ra dX 



CpS N — CpSjY Une Tatm,N{Xline) T N (Xline) Ra (Aline) 

(15) 

If the N band is sufficiently narrow, Ra and T a tm,N can 
be assumed constant within it. Thus: 



EW 



c P s Jv,;iT7 



cps N - cps N line 



f °°T N (X) dX 

Tjv(AHne) 



(16) 



The equivalent width, EW, can therefore be obtained 
from (|12[) and (|16[) without requiring photometric calibra- 
tions. It is equal to the ratio between the count rates of the 



line (cps N line ) and the continuum (cps N — cps^, 



plied by the "N-band effective width" 



SS°T N (\)d\ 



multi- 
The 



Tjv(A H „ e ) 

"N-band effective width" is a measure of the N-band width 
multiplied by the ratio between the N-band peak transmit- 
tance and the N-band transmittance at the emission line's 
wavelength, Trf(\u ne ). 

The EW measurement error contributed by atmo- 
spheric effects originates only from cpsjy;^. In order to 
understand how changes in the atmospheric transmittance 
(e.g., clouds) affect the accuracy of the result, we write ()12|l 
with WNCR in the last term expressed through fundamental 
variables, as defined in ((3j and ([10}: 



°P s JVJi 



cp Sjv 



WNCR/ ' 

T„ t 



T W (A) 



\,cont T atm fy(A^ ne ) Tjv(A Hne ) 



Ra dA 



T„ 



'W T„(A) 



'(Aline) T w (\ Hne ) 



Ra dA 



(17) 

It can be seen here and in (1161) that atmospheric 
conditions affect EW only through the relative transmis- 

T atTOjJV (A) anc 1 T atTT1 W (AJ 



T„ 



/(A, 



e) 



T„ 



'(A; 



») 



"Grey" wavelength- 
independent effects do not change these terms, and therefore 
will not affect the result. 

The effect of clouds blocking the field during parts of 
the exposure will be the same as reducing T a t m ,jv(A) and/or 
T a t m ,vK(A) by an unknown factor. This effect is virtually in- 
dependent of wavelength, and therefore does not contribute 
to the measurement error in EW, whether the W and N 
exposures were taken in similar airmass and atmospheric 
conditions or not. 



3 PHOTOMETRIC CALIBRATIONS 

Unlike EW, measuring the line's flux, Fu ne , requires a pho- 
tometric calibration. Fu ne can be obtained from ((7]), if 
r Fatm,N(Xune) ■ Ra (Mine) is known. A useful method to mea- 
sure T atmi jv(Aiine) -RA(A!i„ e ) utilises photometric calibra- 
tion in the N band. Spectrophotometric standard stars (i.e. 
stars with known Fa) should be measured at several airmass 
(AM) values. The measured count rate of the standard stars 
can be expressed as: 
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p oo 

C P S N spec* {AM) = J F A, spec* T atm,cal 

(\,AM)T N {\)R x d\ 

Jo 

(18) 

where: 

cpSjy apect (AM) is the count rate of the spectrophotomet- 
ric standard star imaged through the narrow-band (N) filter 
as function of the airmass (typically in ADU s^ 1 ); 

Fa, spec* is the spectral flux density of the spectrophoto- 
metric standard star; 

T a tm,cai{^, AM) is the atmospheric transmittance during 
the photometric calibration as function of wavelength and 
airmass. 

If the N band is sufficiently narrow, Ra and T atm ,cal 
can be approximated as constant within the N band. This 
yields: 

T ( \ A M\ R <\ \- CpS N,s P ec {AM) 

J- atm,cal\Aline) AJH J ' (inline) — T°°~F — / \ \ i, 

J r A. spec* ljv(,AJ <ZA 

(19) 

From measurements at different airmasses, a fit can be 
made to the following logarithmic relation: 

-2.5l0g [T a tm.cal(Xune,AM) ■ Rx(Xune)} * AM ■ E Ah „ c + Mo 

(20) 

where: 

Ea Hiic is the extinction coefficient of the atmosphere at 

A; jne i 

Mo is the zero point of the fit. 

If the N band image is taken at similar atmospheric 
conditions as during the calibration, and with AM for which 
the fit ((20)) still holds, T a t m ,N(Mine)-R\(\line) can be simply 
obtained by: 

T„(A !ine ) • R x (X Hne ) = io-«-< AM «- E ^ ln e+ M o) (21) 
where: 

AMjv is the airmass through which the N band image was 
taken. 

When photometric calibrations of the N band cannot be 
obtained, one can calibrate the emission-line's flux based on 
wide-band calibrated images. The target has to be measured 
through several calibrated wide bands, to obtain its flux per 
wavelength, Fa, averaged over each band. From these av- 
eraged Fa values a model for the target's continuum flux 
as function of wavelength, F\ :Cont , can be derived. Using ([6]) 
with the modelled F\ tCont {Xi inf ,), and the measured EW, the 
flux can be simply obtained: 



Fline — EW ■ Fa, cont(Xline) 



(22) 



When available, wide band magnitudes from the Sloan 
Digital Sky Survey (SDSS - Abazajian et al. 2009) can be 
used to model Fa, con*- In some cases a linear fit of the g, 
r and i magnitudes (whose average wavelengths are 4686 A, 
6165A, and 748lA respectively) is a sufficient approximation 
for the continuum at Ha. The error of this linear approxi- 
mation can be estimated using the residuals of the fit. 




6200 6400 6600 6800 7000 

MA) 



Figure 1. FA, CO nt of 1, 100 Myr and 10 Gyr old SSPs at z=0, nor- 
malized so that the average Fa, cont over the R-band of each stellar 
population would be 1. Data adopted from the Bruzual-Charlot 
atlas available in ftp://ftp.stsci.edu/cdbs/grid/bc95/templates/ 
(also used for the calibration of several Hubble Space Telescope in- 
struments). The dotted horizontal line shows the model adopted, 
when calibrating using the R-band alone (F\ cont constant over 
the R band) . The dashed vertical line indicates the He? wavelength 
(Xhc)- 



Since typically the emission-line's wavelength is within 
the W band, one might be tempted to approximate 
F\,cont{Xnne) to be the average continuum flux measured 
for W (when a photometric calibration for band W is avail- 
able). Such an approximation is equivalent to a model in 
which Fa. cont is constant within the W band, and is there- 
fore highly inaccurate as evident from Figure [T] 

Figure [1] shows three examples of normalized Fa, cont 
curves (lMyr, 100 Myr and 10 Gyr old stellar populations). 
The naive model of constant Fa. cont over the R band is also 
shown (dotted horizontal line) . Using a photometric calibra- 
tion for W alone (R in this example) implies that the dashed 
line is used instead of the real value. For example, calibrating 
with R filter photometry alone overestimates Fa, cont (Xhc) 
of the 1 Myr old population at z~0 by 13%. This systematic 
calibration error depends on the age of the stellar population 
as well as on its redshift, and propagates directly through 
to a systematic measurement error in Fu nf .. 



4 THE EFFECT OF ERRORS IN WNCR 

In this section we derive the effect of the WNCR estimation 
error (AWNCR) on the measurement errors of the equiva- 
lent width and flux of the line (AEW and AFi ine ). 

Using (HI} we derive (Acps Nj!ine ) WNCR , the measure- 
ment error of cps N Hne propagated from AWNCR (i.e., 
AcpSjy line if all other sources of error are negligible). For 
reasons discussed above, the right-hand side factor of (|12|l 
typically introduces a small correction to the left-hand side 
factor. For error propagation purposes this correction can 
be ignored. Hence: 



© 2011 RAS, MNRAS 000, [HS] 



Continuum removal in Ha measurements 5 



(Acps JViiine ) 



= AWNCR ■ 
^ AWNCR ■ 



dWNCR 
cps w 



(23) 

WNCR 2 

For the purpose of error analysis, we assume that the 
line's contribution to cps w is minor (i.e. 
This is true in typical conditions, where the W-band width 
is significantly larger than the EW. Using (| 10p with this 
assumption we obtain: 



(Acps Jvi „ le ) wNCR - 
From this and Q, (AF i; 



AWNCR 
WNCR 



. cant 



(24) 



/ncr ' ^ ne measurement er- 
ror in Fii ne propagated from AWNCR, can be simply de- 
rived as: 



(AFii„ e ) WNCR — (Acps A , ;ine ) wNCR ■ 



dFu r 



AWNCR 



CP S JV„ 



dcps N line 



(25) 



WNCR T at 
Using (|13p and the assumption made in (]16[) yields: 

„ AWNCR /,r^(A) rfA 

l^'iMJwSOE- WNCR ' Tiv(A Ilne ) '^^n^inej 

(26) 

(AFiine) WNCR is therefore affected by the relative error 
in WNCR, the N-band effective width, and the underlying 
continuum flux. Note that (AF;ine) WNCR does not depend 
on the line's flux itself (Fune)- 

Similarly, (AEW) WNCR , the measurement error in EW 
propagated from AWNCR, can be derived from (fTr?) 
and ((24)) : 

AWNCR J °°Tn(A) d\ / cps NMne \ 
^ WNCR' T N (X Hne ) A cpB NiCOn J 

(27) 

The relative error, obtained by dividing this equation 
by 1(15). is: 



(AEW), 



(AEW) V 



EW 



AWNCR 
WNCR 



1 + 



CPSJV., 



cpsj, 



(28) 



Since the last term in (|27[) is at least unity, it is evident 
that the relative error of WNCR multiplied by the N-band 
effective width is a lower limit of the uncertainty in EW: 

AEW^.iS^. (29) 



Tjv(A;ine) 



WNCR 



Similarly, it is evident from (|28|) that the relative error 
of WNCR sets a lower limit to the relative error of EW: 



AEW > AWNCR 



EW 



WNCR 



(30) 



5 WNCR ESTIMATION METHODS 

For well-resolved galaxies, the continuum light can be di- 
rectly measured from the galaxy itself. This can be done 



using the pixel-to-pixel method (Boker et al. 1999; Rossa 
& Dettmar 2000; Kennicutt et al. 2008; Finkelman et al. 
2010) where cpsjy is plotted for each pixel as function of 
cps^. To reduce the data scatter one may first reject a user- 
specified percentage of the brightest pixels, suspected to be 
foreground stars. Galaxies with smoothly varying light dis- 
tribution, such as elliptical galaxies, can also be fitted with 
isophotal maps in both the W-band and N-band images. 
From these, one can plot the mean cps N vs. mean cps w for 
pairs of isophotes of the same semidiameter (e.g., Macchetto 
et al. 1996). 

The resulting relation is expected to show most data 
points (presumably lacking line emission) aligning along a 
narrow straight line down to the noise level. The slope of 
this line yields the WNCR estimate. A smaller number of 
data points, which trace the Ha emission, are expected to 
be above the line, whereas a small amount of scattered data 
points are also expected below the line. The latter trace 
either regions of excess continuum emission, or excess line 
absorption. Note that this WNCR estimate is affected by 
widespread diffuse emission and by underlying stellar Ha 
absorption across the galaxy. 

A hidden and problematic assumption in applying this 
method is that of a constant colour. The fitted regions 
may show colour variations caused by changes in stellar 
population across the galaxy and in the line-emission 
regions. These may be related to recent or ongoing SF 
episodes. Note that the derived WNCR estimate, which 
is affected by colour, does not represent the emission line 
region but is rather an average value across the galaxy. 

An alternative method to estimate WNCR is by calcu- 
lating the average cps^, to cps N ratio of foreground stars 
(e.g., Pogge, Owen & Atwood 1992; Rand 1996; Koopman, 
Kenney & Young 2001; James et al. 2004). This method 
allows measuring line emission even if the continuum- 
dominated regions of the target cannot be properly sampled. 

The main assumption is that the foreground stars rep- 
resent the mix of stars in the target galaxy, thus match- 
ing on average the spectral distribution of its continuum. 
In order to represent a continuous distribution function of 
stellar types, as many foreground stars as possible should 
be measured. However, regardless of how numerous, fore- 
ground stars represent the local population of the region 
in the Milky Way through which the observation is made, 
rather than the stellar population in the target galaxy. These 
two populations may be significantly different. 

This method is more useful for target galaxies that 
are sufficiently distant to allow the N-band to transmit 
their rest-frame Ha while blocking the local Ha (6563A) 
of Milky Way foreground stars. If this is not the case, 
applying this method implies also the assumption that 
the fraction of foreground stars with strong Balmer line 
absorption or emission (typically O, B and A stars) is 
insignificant. This assumption typically holds for observa- 
tions in mid to high Galactic latitudes. At these latitudes 
Galactic foreground stars between the observer and a 
target tend to be of types later than A due to the age 
of the halo and the shape of the stellar initial mass function. 

We note that an alternative set of bands may be used, 
in which an off-line narrow band replaces the wide W-band 
(e.g. Almoznino & Brosch 1998; van Zee 2000; Dale et al. 
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2010). If the central wavelengths of the off-line and on- 
line bands are close enough, the WNCR uncertainty due to 
colour effects may be reduced. However, this comes at the 
expense of increased random errors due to the lower depth 
of images obtained with limited telescope time. An alter- 
native approach utilises two medium to wide off-line bands 
flanking the Ha line (e.g. MacKenty et al. 2000; Lin et al. 
2003; Calzetti et al. 2004; Ly et al. 2011). This approach re- 
duces colour effects, while having a smaller impact on image 
depth. 



6 COLOUR EFFECTS 

Colours of the studied emission-line regions may differ from 
those of the foreground stars or those of the continuum- 
dominated galactic regions used for estimating WNCR. As 
a result, the estimated WNCR may deviate from its value in 
the studied region. Furthermore, since colours typically vary 
within a galaxy, different parts of a target galaxy may have 
different WNCR values. In an attempt to improve accuracy, 
the value of WNCR is sometimes adjusted subjectively to 
obtain a desirable continuum subtracted image (Rand 1996; 
James et al. 2004; Kennicutt et al. 2008). 

Although the uncertainties in estimating WNCR are 
typically of few percent, these can result in a much higher 
uncertainty in the measured Ha flux, especially for low EW 
values (Macchetto et al. 1996; Hameed & Devereux 1999; 
Fig. 7 in Koopmann et al. 2001). For example, for a 5% 
error in WNCR and an N-band effective width of 100A, the 
relative error of EW would be at least 5% (eq. 130)1 . while 
the absolute error of EW would be at least 5A (eq. I29p . For 
low SF rates the latter can be comparable to the measured 
EW. 

In order to estimate the variation of WNCR with colour 
we model observations of different objects with recession 
velocities up to 30000km s~ x (z=0.1). For each object we 
simulate observations through a Johnson- Cousins R filter 
and a narrow-band filter, 50A wide and centred on the red- 
shifted Ha line. To simulate foreground stars we use vari- 
ous spectra of Main Sequence (MS) and off-MS stars, taken 
from the Bruzual-Persson-Gunn-Stryker atlas (available at 
ftp:/ /ftp. stsci.edu/cdbs/grid/bpgs/ ). The continua of differ- 
ent parts of target galaxies are represented by redshifted 
spectra of single stellar population (SSP) models at differ- 
ent ages in the range 10 6 -10 10 yr. These star-bursts are char- 
acterised by a Salpeter initial mass function (IMF) with a 
lower and upper mass of 0.1 and 125 Mq, respectively. The 
data are taken from the Bruzual- Chariot atlas (available in 
ftp: / / ftp.stsci.edu/cdbs/grid/bc95 / templates/) . 

Using the spectra of these objects, the filters' transmit- 
tance curves, and the responsivity of a typical electro-optical 
system, we calculate cps^, cont and cps^ cont and from it the 
theoretical WNCR for each object (equations [51 and l 101) . The 
calculated WNCRs for each narrow-band filter are plotted 
vs. the computed SDSS g-r colour of each object (Figure[2]). 
For each narrow-band filter WNCR values are normalized 
using the WNCR of a theoretical object with constant /„ 
(zero colour in the AB magnitude system). Linear fits for 
foreground star points are also shown in each plot. 

As can be seen in Figure [21 WNCR of the foreground 
stars depends almost linearly on g-r. WNCR of the simu- 
lated target galaxy parts seem to follow the curve of the 



foreground stars with the exception of two points (10 7 yr 
old SSP population with recession velocities of 20000 and 
30000 km s" 1 ). 

Based on this, we suggest the following method to ob- 
tain a more accurate estimate of WNCR. A colour measure- 
ment of the objects in the imaged field should be obtained 
either by imaging in a second wide band, or by using public 
databases. WNCR values should be measured for each ref- 
erence object (foreground star or part of the target galaxy 
with no line flux) by dividing the counts in W by the counts 
in N (equation I10p . A linear relation should then be fitted 
to the measured (colour,WNCR) points. Finally, this linear 
relation should be used to obtain a WNCR estimate for each 
part of interest in the target galaxy as function of its colour. 

Let us now estimate the possible errors in WNCR, 
when obtained without colour information. First, we esti- 
mate the normalized WNCR value obtained using the fore- 
ground stars method. For this we measure the colour of fore- 
ground stars around arbitrary nearby galaxies. This is done 
by collecting data from the SDSS archive for regions around 
these galaxies, each with a radius of 30 arcmin, and con- 
taining ~10,000 objects classified as stars. The g-r colour 
histograms for two such regions are presented in Figure 
[3] Both histograms show a distribution that peaks at g-r 
of ~0.8, corresponding to lower main sequence and giant 
branch stars. Since the average colour is also of a similar 
value, we adopt 0.8 as a typical g-r colour for foreground 
stars around galaxies. 

The error in WNCR can be estimated for the different 
recession velocities of Figure [5] by comparing the WNCR 
value of foreground stars at g-r = 0.8 with the WNCR val- 
ues of the simulated target galaxy parts (the redshifted stel- 
lar populations). When, for example, measuring a starburst 
galaxy dominated by a young stellar population (10 6 yr) 
with foreground stars of g-r = 0.8, WNCR will be under- 
estimated by up to ~30% depending on the galaxy's red- 
shift (for cz 30000 km s" 1 ). This means underestimating 
the Ha EW by a similar fraction (equation I30[) . and under- 
estimating Fii ne by ~0.3 of the continuum flux within the 
N-band (equation I26|) . For older underlying stellar popula- 
tions the error in WNCR would be smaller. For 10 9 yr the 
underestimation will reach 8% at 10000 km s _1 , while for 
10 10 yr WNCR will be overestimated by up to ~10%. 

Estimating WNCR using galaxy parts believed not to 
form stars may also be problematic, since it does not take 
into account the possible colour variations across the galaxy 
and the expected colour change due to star formation in the 
Ha emitting regions. Let us consider a case where a 10 6 yr 
old stellar population is observed within an elliptical galaxy 
dominated by a 10 10 yr old population. The flux from the 
young population is expected to overwhelm the flux from 
the old population, even if it is only a small fraction of the 
entire star population. Therefore, the measured WNCR will 
correspond to redder objects, i.e. would be underestimated 
by up to ~40% depending on redshift (for cz ^ 30000 km 

These errors in WNCR propagate through Fu ne to SF 
measurements. Since the errors depend on colour properties 
of the galaxies, a bias in the colour properties of a sample 
may result in a bias in the SF properties concluded by a 
survey. 

This applies also to surveys that measure cosmic SF 
rate at low redshifts (low z points in the Madau plot). Since 
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Figure 2. The theoretical WNCRs measured for a set of objects with a set of narrow-band filters. Foreground stars of various spectral 
types and luminosity classes are represented by X , while redshifted SSPs evolving with time are marked by filled shapes ( • - 10 10 yr, 
■ - 10 9 yr, T - 10 8 yr, ♦ - 10 7 yr, A - 10 6 yr). The set of filters sample the rest-frame Ha line in galaxies with recession velocities of 
5000, fOOOO, 20000 and 30000 km s" 1 . The WNCR values measured with each narrow-band filter are normalized by the corresponding 
WNCRo value calculated for a constant /„ spectrum (zero colour in the AB magnitude system). The solid lines show linear fits of the 
foreground stars. 
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Figure 3. The distribution of the g-r colour of stars in two regions in the sky. The regions were arbitrarily selected to include all stars 
in a circular aperture of 30arcmin around two arbitrary nearby galaxies (NGC 5363 and NGC 7625). 
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star forming galaxies tend to be bluer than foreground stars, 
results may be biased. As discussed above, such a survey that 
uses foreground stars and no colour information may result 
in the underestimation of WNCR by tens of percents. Such 
underestimation, which increases with redshift, will result in 
the underestimation of the cosmic SF rate. 



7 CONCLUSIONS 

We identified an important error in estimating extragalactic 
Ha emission, originating from commonly used procedures 
of continuum subtraction. We showed that when measuring 
an emission line by observing through a narrow-band (N) 
and a wide-band filter (W), where the N band contains the 
emission-line, and the W band may or may not contain it: 

(i) The Equivalent Width, EW, can be obtained without 
the need for photometric calibrations using (1121) and (|16[) . 
To improve accuracy: 

• the N filter should have a very narrow band and neg- 
ligible "out of band leaks" ; 

• the exposures of the W and N bands must be taken 
in very similar atmospheric conditions and air-masses, or 
with measured atmospheric transmittance at the line's 
wavelength [T a tm,w{^iine) and T a t 

(ii) An accurate measurement of the line flux, Fjjne, can 
be obtained by photometric calibrations in the N-band, if it 
is sufficiently narrow to assume that the responsivity, Ra, 
and the atmospheric transmittance, T a t m ,jv, are constant 
within it. Alternatively, if calibrated measurements of 
the magnitude in several wide-bands are obtained, a fit 
can be used to estimate the continuum flux under the 
line, Fx Une ,cont, from which Fu ne can be obtained. Using 
calibrations in the W-band alone may lead to systematic 
errors of over 10%. 

(iii) An accurate estimation of WNCR, the W to N band 
ratio of count rates contributed by the continuum (equation 
I10p . is crucial for accurate emission line measurements: 

• The uncertainty in EW is equal to or greater than the 
relative error in WNCR multiplied by the N-band effective 
width (equation I29|) : 

• The relative uncertainty in EW is equal to or greater 
than the relative error in WNCR (equation 1301) : 

• The uncertainty in the emission- line flux, Fu ne , is 
equal to or greater than the relative error in WNCR mul- 
tiplied by the underlying continuum flux and by the N- 
band effective width (equation I26|) . 

(iv) The value of WNCR may vary significantly with 
the colour of the measured object. For typical Ha mea- 
surements at redshifts cz ^ 30000 km s -1 the dependence 
on g-r was found to be almost linear. Where foreground 
Milky Way stars are used for estimating WNCR, colour 
effects may lead to an underestimation as great as 30%, or 
in other cases overestimation by as much as 10%. Where 
regions in the target galaxy itself are used for the esti- 
mation, WNCR may be underestimated by as much as 40%. 

(v) We suggest estimating WNCR as a linear function 
of target colour, rather than as a constant value. This will 
reduce errors, and will make the measurement procedure 



more objective. 

(vi) The errors described above can bias the results of a 
study of a class of galaxies by systematically underestimat- 
ing or overestimating the emission-line flux. Cosmic SF rate 
measurements at low redshifts may also be biased by this 
effect. If not corrected properly, such measurements may 
underestimate the low z points in the Madau plot. 
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